Backreaction in Acoustic Black Holes 
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The backreaction equations for the linearized quantum fluctuations in an acoustic black hole 
are given. The solution near the horizon, obtained within a dimensional reduction, indicates that 
acoustic black holes, unlike Schwarzschild ones, get cooler as they radiate phonons. They show 
remarkable analogies with near-extremal Reissner-Nordstrom black holes. 
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One of the most surprising and far reaching result for 
its implications in modern theoretical physics is the pre- 
diction made by Hawking [1] that black holes emit ther- 
mal radiation at a temperature Th proportional to the 
surface gravity k of the horizon. For a Schwarzschild 
black hole of mass M, k = (4M)" 1 and T H = ^(SttM)" 1 
(we have set the velocity of light and Boltzman constant 
equal to one). Hawking obtained this result using quan- 
tum field theory in curved space, a scheme for dealing 
with the matter-gravity system where matter is quan- 
tized according to quantum field theory whereas grav- 
ity is treated classically according to Einstein General 
Relativity. The scale at which this framework becomes 
unreliable is the Planck length where the description of 
spacetime as a continuous differcntiable manifold proba- 
bly breaks down. Coming back to black holes, because 
of the quantum emission, they are unstable. Extrapo- 
lating Hawking's result (which is strictly valid only for 
stationary or static black holes) one can conjecture that 
as the mass decreases, the hole gets hotter and hotter (be- 
ing the temperature inversely proportional to the mass) 
and eventually disappears in a time scale of the order of 
the initial mass to the third power. A more quantitative 
analysis can be performed by looking at the first order (in 

K) corrections g^l to a classical black hole metric g^l in- 
duced by the quantum emission. These can be calculated 
using the semiclassical Einstein equations [2] 
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Here is the Einstein tensor evaluated for the quan- 
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the perturbation g£l (of order H). The r.h.s. represents 
the expectation value of the stress tensor for the quan- 
tum matter field evaluated in the classical background 
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In a very interesting paper, appeared in 1981, Unruh [3] 
showed that a thermal radiation similar to the one pre- 
dicted by Hawking for black holes is expected in a com- 
pletely (at first sight) different physical scenario, namely 
a fluid undergoing hypersonic motion. This opened the 
way for the study of condensed matter analogues of 
Hawking radiation [4] , a rather promising field of research 
where the connections to the experimental side do not 
seem so remote, compared to gravity. 

The Eulerian equations of motion for an irrotational 
and homcntropic fluid flow can be derived from the action 
[5], [6] 



S = - 
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where p is the mass density, ip the velocity potential, i.e. 

u the internal energy density and a dot means 
time derivative. 

Varying the action with respect to ip one gets the con- 
tinuity equation 



p + ^.(pV) = 0, (3) 
whereas variation with respect to p gives Bernoulli's 

du > 
dp> 

— >2 
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equation {p, = ^) 



Expanding the fields around a solution (po, Vo) of the 
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FIG. I: A Laval nozzle. The waist of the nozzle represents 
the sonic horizon (|wo| = c). In the region on the right of the 
waist \vq\ < c and on the left \vo \ > c (sonic hole). 



classical equations of motion (3,4), Unruh showed that 
the quadratic action 52 for the fluctuations of the ve- 
locity potential can be written in a simple and elegant 
geometrical form, namely 
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where ipi is the fluctuation and gffi is the so called acous- 
tic metric 
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expressed in terms of the background quantities, c is 
the sound speed, i.e. c 2 = pdp/dp, and / is the three- 
dimensional identity matrix. 

As it can be seen, S2 has exacly the same form of an ac- 
tion for a masslcss, minimally coupled, scalar field prop- 
agating in a "curved spacetime" whose line element is 
ds 2 = g^}dx^dx v . The region of the fluid for which 
v 2 > c 2 is called sonic black hole: its boundary \vq\ = c 
defines the sonic horizon. Sound waves cannot excape 
from this region, since they are dragged by the fluid. 
A typical example is the Laval nozzle of Fig.l. The fluid 
flows from right to left. At the waist of the nozzle the 
fluid velocity reaches the speed of sound: this is the lo- 
cation of the sonic horizon for free fluid motion. 

Using Hawking's arguments Unruh, quantizing the 
field tpi, showed that in the formation of a sonic hole 
one expects a thermal emission of phonons at a temper- 
ature Tjj = hk/(2irc), where k is the surface gravity of 
the sonic horizon [7] 
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n is the normal to the horizon. 

In this paper we will give a first qualitative analysis 
of the effects this emitted radiation has on the fluid dy- 
namics, i.e., the backreaction of the linearized phonons 
in sonic black holes. 

Using standard background field formalism we write the 



fundamental quantum fields p, tp as the sum of back- 
ground fields p, tp (not necessarily satifying the classi- 
cal equations of motion), plus quantum fluctuations, i.e. 
p = p + 8p, tp = tp + Sip. Integrating out the quantum 
fluctuations one obtains the one loop effective action for- 
mally defined as T = S + \h trln[D s ( p ^)] +0(H 2 ), where 
S is the action (2) and O g ( p ^ is the covariant D' Alam- 
bertian calculated from an acoustic metric g^ v {p,v) of 
the functional form of eq.(6) with po and vo replaced by 
p and v respectively [8]. We assume that divergences in 
the determinant of the above effective action are removed 
employing a covariant regularization scheme. This is the 
key hypothesis of our work. We will comment on it later. 
Therefore one can write T = S + S q (g t _ iu (p, v)), where the 
quantum part of the action T depends on the dynamical 
variables p, v via the acoustic metric only, and coincides 
with the effective action for a massless scalar field propa- 
gating in a spacetime whose metric is g^i^p, v). Using the 
chain rule we write ^ L = and similarly for ^S- . 

We write p = p + p q and v = v + v q , where p , v satisfy 
the classical equations of motion (3,4) and p q ,v q are the 
corrections of order O(h) induced by the quadratic fluc- 
tuations of the phonons field. The linearized backreac- 
tion equations, analogues of semiclassical Einstein eq.(l), 
read (assuming for simplicity sake a constant velocity of 
sound c) 



Pq + ^i(PqVOi) + Vj (pn{v qi ~ 

4>q + V ■ Vg + ^ Pq 
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where 



2 SS q 



is the expectation value 



of the quantum version of the so called "pseudo en- 
ergy momentum tensor" T^ v = —-jL=jj£ [6] and (T) = 

gWv" (T^v) . From the above construction it follows that 
{T^iP) coincides with the expectation values of the stress 
tensor for a massless scalar field propagating in an ef- 
fective spacetime whose metric is gffi (p , v ). These ex- 
pectation values should be taken in the sonic analogue 
of the Unruh state [9] (the quantum state appropriate to 
describe black hole evaporation at late times), in which, 
in the remote past prior to the time dependent forma- 
tion of the sonic hole, the quantum field is in its vacuum 
state. Inspection of the backreaction equations reveals 
that eq.(8) is the first order in H conservation equation for 
the Noether current associated to the ip — > ip+ const sym- 
metry of the effective action T. Note that the phonons 
contribution in the Noether flux coming from S q {g^ v ) is 

just the pseudo momentum density \fg~W{Tf). Eq.(9) is 
the first order Bernoulli equation modified by the pres- 
ence of the trace term which represents an additional 
contribution to the chemical potential induced by the 
fluctuations which, as we shall see, causes the fluid to 
slow down. These equations reflect the underlying two- 
component structure of the system as in the Landau- 
Khalatnikov theory of superfluidity [10]. One can also 
rewrite the equation (8) as the usual fluid continuity 
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equation by a simple redefinition of the velocity field 
Vqi — > v q i — (^jr- + -^-{Tjii^j which incorporates the 
phonon momentum density. Consequently the Bernoulli 
equation (9) rewritten in terms of this redefined velocity 
contains, besides the trace, also other terms. 
Unfortunately no explicit solutions of eqs.(8,9) can be 
given since {T^ v ) is unknown. 

In the black hole case, where similar difficulties arise, 
a qualitative insight in the evaporation process can be 
gained using 2D dimensional models [11]. The most pop- 
ular is the one proposed by Callan, Giddings, Harvey and 
Strominger (CGHS) [12] in which the 4D quantum stress 
tensor is replaced by a 2D one associated to a minimally 
coupled massless scalar field described at the quantum 
level by the Polyakov action [13]. 

With the same spirit a qualitative description of the back- 
reaction in a hypersonic fluid can be obtained assuming 
a one dimensional flow for the fluid, let's say along the 
axis of the nozzle, the z direction. So all physical fields 
will depend only on t and z. 

The effective action for the CGHS-likc model for the fluid 
quantum dynamics can be then given as 



r< 2 ) = + s, 
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is obtained integrating S over the transverse coordinates 
x, y and A is the area of the transverse section of the 
nozzle. S po i is the Polyakov action. The backreaction 
equations following from r' 2 -* are 



Apq 



d z [A(p q v + p Vg)] 
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(T^ b ') is the quantum stress tensor for a massless scalar 

field minimally coupled to the (t, z) section of the 

acoustic metric gffi of cq. (6). The trace (T^) is com- 
pletely anomalous and given by 

(T(2)> = i E(2) ' (14) 

where is the Ricci scalar for the metric . The 
phonons expectation values appearing in the conserva- 
tion equation (12) can be easily expressed by a coordinate 
transformation in terms of (T±±) where = t±z± and 
z* ± = J dz/(c =F t>o) -1 ) f° r which the Polyakov approxi- 
mation gives 



(t±±) = -h{i2ir)- l c 1,2 c;±i; + a 



± • 



(15) 



Here C is the conformal factor for the 2D metric, C — 
Po(c 2 — Vq)c~ 3 , and for the Unruh state A + = 0, A_ = 



hk 2 /48irc 4 . 

Assuming as profile for the Laval nozzle A = A 
with A, (3 constant, the classical solution reads 



Po = 



pe 
A 



~c? , p = const 
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where the sound velocity c is taken to be constant. The 
location of the sonic horizon is z = 0, where v = — c, and 
the region z < is the sonic black hole for which \vo\ > c. 
These expressions should be regarded as describing the 
classical would be asymptotic configuration of the fluid 
resulting from the (time dependent) formation of a sonic 
hole [14]. 

Expanding the stress tensor and the background quanti- 
ties near the horizon z = one eventually arrives [15] to 
the following solution for the velocity 

V = V() + V q ~ — C + CKZ — ^K 2 Z 2 + e(6l + C\KZ)Kt , (17) 

where e = ft/(j4 2 pe _1//2 c) is the dimensionless expan- 
sion parameter and b\ — 9-f/2, c\ = —3047/15, where 
7 = Ac 2 k 2 /24ir, and k — \J (3 / A has dimension L _1 and 
is related to the surface gravity k — k/c 2 . 
The solution is valid for kz -C 1 and cut 1. 
The boundary conditions imposed on the backreaction 
equations have been chosen so that at some given time 
(say t = 0) the evaporation is switched on starting from 
the classical configuration (u ,/3o), i-e. p q {z,t = 0) = 
v q {z,t = 0) = %(z,t = 0) = 0. 

Inspection of eq.(17) shows the net effect of the backre- 
action. Being kz <§C 1, from eq.(17) we have v q > 0, i.e. 
the fluid is slowing down. This goes with an decrease of 
the density (p q < 0) in the same limit [15]. 
Eq.(17) allows also to follow the evolution of the acoustic 
horizon. Being the horizon defined by v = — c, this yields 



ZH ^ 



ebit 



(18) 



i.e. the horizon is moving to the left with respect to the 
nozzle: the hypersonic region shrinks in size. The coef- 
ficient b\ determining the quantum correction to the ve- 
locity (17) and hence the evolution of the horizon is just 
the gradient of the additional chemical potential related 
to the expectation value of the trace evaluated at z = 0. 
This should be compared to the black hole case where 
the evolution of the horizon is determined by the energy 
flux (M oc (T r t ) in spherical symmetry [2]). While in 
the latter case Hawking radiation occurs at the expense 
of the gravitational energy of the black hole, in the fluid 
phonons emission takes away kinetic energy from the sys- 
tem. 

From cqs.(17, 18) one can evaluate the correction to the 
emission temperature 



h dv 
Tu= 2^dz 



he 

27' 
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Using particular values for liquid helium, the fractional 
change of the temperature per unit time is of order 
10- 8 a- 1 . 

The expression we have obtained is rather significative. 
Unlike a Schwarzschild black hole the emission temper- 
ature for a sonic black hole decreases in time as the ra- 
diation emission proceeds. This behaviour is reminis- 
cent ^of the near-extremal Reissner-Nordstrom black hole 
(M > |Q|) where Q is the conserved charge of the hole. 
As the mass decreases because of Hawking evaporation, 
the Hawking temperature decreases as well, vanishing 
when M = \Q\. This ground state is approached in an 
infinite time (third law of black hole thermodynamics). 
One can conjecture a similar behaviour for the sonic hole 
with a vanishing temperature ground state. 
The basic question which remains open is to what extent 
our results do depend on the dispersion relation assumed 
(free scalar field), which ignores short distance correc- 
tions due to the molecular structure of the fluid, and on 
the covariant regularization scheme used to subtract ul- 
traviolet divergences. These two aspects are deeply con- 
nected. For the Polyakov theory we have used in the 2D 
backreaction equations (12, 13), Jacobson [16] has argued 
that, within a covariant regularization scheme, no signif- 
icant deviation from the usual expression for the trace 

( 2) 

anomaly and the components of {T ab ) are expected if 
one introduces a cutoff at high frequencies. However, for 
the hydrodynamical system we have considered, covari- 
ance is a symmetry of the phonons low energy effective 
theory only, which is broken at short distance. Hence non 
covariant terms depending on the microscopic physics are 
expected to show up in the effective action and are crucial 
for a correct description of the unperturbed quantum vac- 
uum of the fluid. However the expectation values (T^) 
entering the backreaction equations (8, 9) do not repre- 



sent the energy momentum of the fluid quantum vacuum. 
They describe instead the perturbation of the stationary 
vacuum (whose energy is strictly zero [17]) induced by 
inhomogcncities and by the time dependent formation of 
the sonic hole which triggers the phonons emission. In 
this paper we have assumed that these deviations can be 
computed within the low energy theory. This situation 
is not unusual. Casimir effects are well known exam- 
ples of vacuum disturbances caused by the presence of 
boundaries. It happens that the Casimir energy is often 
(but not always, see G. Volovik in Ref. [4]) independent 
on the microscopic physics and can be calculated within 
the framework of the low energy theory. This happens 
because, while low frequency modes are reflected by the 
boundaries, for the high energy ones the wall is transpar- 
ent. They produce a divergent contribution to the vac- 
uum energy which is canceled by a proper regularization 
scheme and does not affect the finite result. We have as- 
sumed that a similar decoupling happens for the acoustic 
black hole. The check of our hypothesis would require an 
analysis of the quantum system within the microscopic 
theory which takes into account the time dependent non 
homogeneous formation of the sonic hole. This is for the 
moment beyond computational capability. Anyway, it 
has been shown that modifications of the dispersion re- 
lations, to take into account short-distance behaviour of 
the high-energy modes, basically do not affect the spec- 
trum of the emitted phonons [18]. This is not a proof 
that observables like (T^ v ) are also unaffected by short- 
distance physics. However, it can be an illuminating hint 
taking in mind what (T^„) does really represent and the 
indications coming from the Casimir effect. 
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